Jamming under tension in polymer crazes 
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Molecular dynamics simulations are used to study a unique expanded jammed state. Tension 
transforms many glassy polymers from a dense glass to a network of fibrils and voids called a craze. 
Entanglements between polymers and interchain friction jam the system after a fixed increase in 
volume. As in dense jammed systems, the distribution of forces is exponential, but they are tensile 
rather than compressive. The broad distribution of forces has important implications for fibril 
breakdown and the ultimate strength of crazes. 

PACS numbers: 61.43.Fs, 62.20.Fe, 83.10.Rs 



The common features of jammed systems ranging from 
molecular glasses to granular media have sparked great 
interest [1,2]. These systems jam as the available volume 
becomes too small to allow relative motion of particles. 
Unifying features of the jammed state are the presence 
of an exponential force tail at large compressive forces 
and the appearance of a network of forces on scales much 
larger than the size of the constituent molecules or grains 
[3,4]. The origins of these features are still debated. 

In this Letter, we consider a qualitatively different 
jammed state that forms under tension. Many amor- 
phous polymers expand from a typical dense glass to a 
"craze" consisting of a network of fibrils and voids [5] . As 
shown in Fig. 1, the dense and expanded jammed states 
coexist at a fixed tensile stress S, much like liquid and 
gas phases coexist at a fixed pressure. There are no cova- 
lent bonds that keep the craze from unravelling, but ex- 
periments suggest that the topological constraints called 
entanglements that limit dynamics in polymer melts be- 
have like crosslinking bonds [5]. These experiments can- 
not address how deformation occurs, the distribution of 
forces within the system, or the configurations of indi- 
vidual chains within the intricate craze structure. We 
describe extensive molecular dynamics simulations that 
address all of these issues. We find that the distribution 
of tensile forces in the craze has an exponential tail anal- 
ogous to that observed for compressive forces in dense 
jammed systems. The polymer undergoes an approxi- 
mately affine displacement as it deforms into the craze. 
Expansion stops when segments that are only 1/3 of the 
distance between entanglements are pulled taut. This 
factor can be understood from simple geometric argu- 
ments and the assumption that entanglements act like 
chemical crosslinks between chains. 

We study a standard coarse-grained model [6], where 
each linear polymer is modeled by N beads of mass 
to. Van der Waals interactions are modeled with a 
standard 6-12 Lennard- Jones (LJ) potential with energy 
scale e and length scale a. A simple analytic poten- 
tial [7,8] is used for covalent bonds, with equilibrium 
length Iq = 0.96fi between adjacent beads along the 



chain. In this fully flexible (fl) model, the number of 
beads between entanglements (entanglement length) is 
w 70 [6]. In order to analyze entanglement effects, 
we also consider semiflexible chains (sfl) where an addi- 
tional bond-bending potential stiffens the chain locally 
and reduces the entanglement length to Nf 1 ss 30 beads 
[8-10]. We consider two temperatures T = O.le/fcs and 
T = Q.3e/kB, where the latter is close to the glass transi- 
tion temperature. The strength of adhesive interactions 
between beads is varied by truncating the LJ force at 
either r c = 1.5a or r c = 2.2a [11]. 




FIG. 1. Close-ups of interface between dense polymer and 
craze for (a) T = 0.1e/k B , r c = 1.5a (fl) and (b) T = 0.3e/fc s , 
r c — 2.2a (sfl). The lateral dimension of both panels is 64cr. 
The diameter and spacing of the fibrils increase as T and/or 
r c increase. 



An amorphous glassy state is created in a cubic simu- 
lation cell of side L using standard techniques [6]. The 
period in the z-direction is then increased at a small con- 
stant velocity [11]. After an initial nucleation period, 
the system separates into two phases (Fig. 1). The craze 
network grows at constant tensile stress S through plas- 
tic flow in a narrow interfacial region called the "active 
zone." The volume occupied by the polymer increases 
by an "extension ratio" A that is remarkably insensitive 
to all parameters other than N e . For example, increas- 
ing the temperature and interaction range leads to much 
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larger fibrils in Fig. 1(b) than (a), yet there is no measur- 
able change in A. From the ratios of the initial and final 
densities we obtain Ag = 6.0 ± 0.6 and A s a = 3.5 ± 0.3, 
independent of N, T, and r c . Experimental values of A 
range from about 2 to 7 [5] . 

Although strain rate is strongly localized at the in- 
stantaneous position of the active zone, the net effect is 
a nearly uniform or affine expansion of each region of the 
polymer. As shown in Figure 2(a), the final height z/ 
of beads with an initial height of Zi is very close to Xzi. 
The rms variation in the final height (half the errorbar 
length), has a relatively small constant value. For rea- 
sons discussed below, it is close to A e /3 for both flexible 
and semiflcxible chains at all T and r c studied. There are 
also much smaller rms displacements in the x — y plane. 
These allow chains to lower their free energy by gath- 
ering into fibrils whose local density is close to that of 
the dense glass. The magnitude of lateral displacements 
varies with the size and spacing of fibrils, which depends 
on both T and r c . For example, the lateral displacement 
increases from <~ 2.5a in Fig. 1(a) to <~ 5.6<r in Fig. 1(b). 

The affine nature of the deformation along z can also 
be inferred from the change in the conformation of in- 
dividual polymers. In the initial state, chains exhibit 
an ideal random walk (RW) structure inherited from the 
melt. The rms distance Ar(AN) between beads that 
are AN neighbors apart scales as Ar 2 = I p IqAN, which 
defines the persistence length l p . In the dense glass, 
the mean-squared projection along each of the three 
axes is equal. An affine deformation by A along z will 
only change the displacements along the z-axis, yield- 
ing (Az 2 ) = X 2 l p l Q AN/3. Fig. 2(b) shows that the final 
conformation of chains follows this expression at large 
scales. The in-plane components of the displacement are 
not changed during crazing. 




zjo- AN 
FIG. 2. (a) Final vs. initial height for flexible (large 
slope) and semiflexible (small slope) chains (T = O.le/fcs, 
r c — 1.5a). Averages were performed over z-intervals of width 
a. Straight lines have slope Aa = 5.9 and A s a = 3.5, respec- 
tively. Error bars show the standard deviation from the av- 
erages at each height and are 17 ± la (fl) and 9 ± la (sfl). 
(b) Height change Az as a function of distance AN from the 
chain center. Straight lines have slope A 2 / p Zo/3 with A from 
(a). Systems at different T,r c and N show the same results. 

The expansion along z pulls short segments of the 



chains taut, so that Az 2 rises quadratically at small 
A A in Fig. 2(b). The typical number of beads in a 
taut segment, A st , can be obtained from the intersection 
between the two asymptotic scaling forms: (N st lo) 2 = 
X 2 l p loN st /3, yielding A st = X 2 l p /3lo- Inserting the ob- 
served values of A and l p (l p = 1.65a and l p R = 2.7c), we 
arrive at A s fl t ~ 21 ± 4 and A s f ~ 12 ± 2, respectively. 
The length of taut sections can also be determined by 
direct analysis of the chain geometry. Fig. 3(a) shows 
the distribution P(A st ) of the number N st of successive 
beads whose displacement continues upwards or down- 
wards. Here a bond is considered up (down) if it is within 
45° of the +z (-z) axis. For both flexible and semiflexible 
chains, P(N st ) has an exponential tail with a characteris- 
tic decay length that, like A, is independent of A, T, and 
r c . The decay lengths, A s fl t ~ 24±3 and A s f ~ 13±2, are 
in good agreement with the prediction from RW statis- 
tics. Essentially the same length scales appear in the 
decay of the correlation function for the z-component of 
successive bonds, Fig. 3(b), or the bond-angle correlation 
function (not shown). 

It is not surprising that all the above lengths are com- 
parable to the deviations from a purely affine deformation 
in Fig. 2(a), since segments of length A st are stretched 
taut from their initial RW configuration. However, a very 
successful expression for A is usually derived by assum- 
ing that the taut segments contain A e beads rather than 
A e /3 [5]. The maximum extension ratio, A max , is defined 
as the ratio between the fully stretched length N e lo and 
the initial end-to-end distance of a RW of N e steps 

A max = NJo/ilolpNe) 1 / 2 = ,JloN e /l p . (1) 

Calculated values of A max are very close to values of A 
measured in experiments and in our simulations where 
Amax = 6.5 ± 0.5 and 3.5 ± 0.3 for flexible and semiflcx- 
ible chains, respectively. However, it was noted in early 
work that fully stretched chains would actually yield an 
even larger extension ratio because the initial end-to-end 
vector of segments of length A e is randomly oriented 
[12]. Since the mean-squared projection along any di- 
rection is only 1/3 of the total, A would be \/3A max for 
fully stretched chains. The observation that A w A max 
implies that the average length of stretched segments is 
only N e /3, as we have shown. However those segments 
that are initially along the z-axis are fully stretched when 
A = A max , and it appears that these few segments are suf- 
ficient to prevent further elongation. Thus the factor of 
1/3 results from the random nature of the entangled net- 
work. Note that a recent study of chains with random 
crosslinks [13] is consistent with Eq. (1), while all seg- 
ments between crosslinks were fully stretched in a similar 
study of ordered networks [14]. 
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FIG. 3. (a) Probability of straight segments of length N st 
for flexible and semiflexible chains. Thick lines: T = 0.1e/fcs, 
r c = 1.5a, N = 512 and 1048576 beads. Dotted lines: 
T = 0.3e/fc B , r c = 2.2a, N = 512 and 262144 beads. Long 
dashed lines: T = 0.1e/k B , r c = 1.5a, N = 256 and 262144 
beads, (b) correlation function for the z-components of the 
bonds for the same systems. Thin solid lines in both pan- 
els show fits to an exponential decay with indicated decay 
lengths. 

A crosslinked system will jam when there is a network 
of fully stretched covalent bonds, but it is less clear how 
entanglements can lead to a jammed network. The entan- 
glement length is typically determined from the response 
of a polymer melt to a sudden strain [6,15]. The resulting 
time-dependent stress shows a broad plateau where the 
polymer has relaxed on short length scales, but is unable 
to relax at larger length scales because interpenetrating 
loops of polymer cannot pass through each other. The 
entanglement length corresponds to the typical length 
along the backbone between these constraints. The stress 
ultimately decays by the slow diffusion of polymers along 
their length until their ends pass through the loops and 
release all remaining constraints. Expansion of the poly- 
mer into a craze is also limited by the inability of polymer 
loops to pass through each other and interchain friction 
prevents the chains from unraveling by the snake- like mo- 
tion that allows diffusion in the melt. The tension pulling 
two interpenetrating loops in opposite directions creates 
large compressive forces at the entanglements. As in con- 
ventional jammed systems this increases the barrier for 
sliding. Long chains have many interpenetrating loops, 



and the tensions pulling in each direction nearly cancel. 
Chains with N < 2N e disentangle rather than forming a 
craze [11]. 

The distance between entanglements is determined by 
the chain statistics [16] which are inherited from the melt. 
They are not sensitive to temperature because there is no 
evolution of large scale structure in the glassy phase, and 
they are relatively insensitive to the strength of adhesive 
interactions and density. Thus A is nearly constant, while 
the fibril structure at scales less than N e may vary dra- 
matically (Fig. 1). We have verified that A only depends 
on chain statistics by artificially varying the persistence 
length l p in the glassy state. As expected, the value of 
A depends only on l p and not on the flexibility of the 
chains, r c , or other details of the interaction potential. 

As in other jammed systems, the distribution of forces 
in the craze follows a universal curve with an exponen- 
tial tail (Fig. 4). The tensions / in the covalent bonds 
along the chains carry most of the stress. Only 10-20% 
of covalent bonds are under compression (/ < 0), and 
this part of the distribution does not change much dur- 
ing crazing. The tensile (/ > 0) part of the distribution 
is exponential, and results for all systems collapse after 
normalizing by the average (/) over the tensile region. 
Intuitively, one might expect the crazes with the larger 
extension ratio (and N e ) to exhibit a larger average ten- 
sion than the less stretched crazes. However, we find that 
(/) is a function of temperature T and cutoff range r c 
only (see caption) and does not depend on N e . Expand- 
ing the craze elastically beyond A increases (/), but the 
normalized probability is unchanged until bonds begin to 
break. 
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FIG. 4. Distribution of tension / in the craze for the sys- 
tems of Fig. 3. The straight line is exp (—//(/}). Only 
bonds under positive tension were included in calculating (/) . 
(/) « 8.2e/a for T = 0.3e/fc B , r c = 2.2a, (/) w 7.2e/a in all 
other cases. 

In equilibrium, the energy of any region follows the ex- 
ponential Boltzmann distribution because this maximizes 
the number of microstates available to the entire system 
at fixed total energy. In jammed systems the total force is 
conserved, and one may imagine that the exponential dis- 
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tribution of local forces also arises because it maximizes 
the number of microstates. Two more detailed models 
for the exponential distribution have been suggested for 
dense jammed systems. 

The "q- model" [4] and its generalizations [17] assume 
that stress propagates unevenly through the system. The 
fraction of stress passed to each neighbor, q, is chosen at 
random from a distribution R(q). For a wide range of R 
one obtains an exponential force distribution. In addi- 
tion, the stress is concentrated in force chains like those 
seen in experiments on granular systems. The fibril net- 
work in a craze naturally provides a branching path for 
transmission of stress. Force is conserved along each fibril 
and then redistributed at the nodes where fibrils merge or 
split. Study of the stress along chains shows that tension 
is correlated along straight segments and then changes 
when the chain changes direction at a node. A more 
detailed comparison to the q-model is in progress. 

Recently, O'Hern et al. [2] have proposed an al- 
ternative explanation for the exponential force tail in 
jammed materials. In equilibrium, the separation of 
beads at small distances is weighted by a Boltzmann- 
factor exp[— V(r)/fcsT], where V(r) is the interaction 
potential. The corresponding distribution of forces is also 
nearly exponential if the force varies rapidly with r. The 
repulsive part of the LJ potential satisfies this criterion, 
but our covalcnt bond potential does not. We have cal- 
culated the distribution of bond energies in the craze and 
find that it is not exponential. 

The presence of an exponential stress distribution has 
consequences for the ultimate fracture of craze fibrils. In 
a minimal model often employed in the literature [18], 
the maximum breaking force -F max of a fibril composed 
of n strands is estimated by assuming that each strand 
carries an equal share of the load and all covalent bonds 
break at a critical force f c . This implies F max / f c = n, 
and all strands break at the same time. If the distribution 
of forces among strands is exponential, the most stressed 
strand will break at a lower F max . The force needed to 
initiate failure can be estimated using a simple scaling ar- 
gument. The first strand will break when nP(f > f c ) = 
n $Tc V(/) ex P(-//(/)) d / = «exp(-n/ c /.F max ) = 1, 
where we have used (/) = F miix /n. This implies a break- 
ing force of F max // C = nj ln(n) instead of n. The re- 
distribution of load after a strand breaks will cause the 
remaining strands to break sequentially. Note that no 
chain scission occurs in the simulations presented here, 
but it does occur upon further straining of the craze. 
A detailed analysis of fibril breakdown will be presented 
elsewhere. 

In summary, we have shown that crazing transforms a 
conventional dense jammed state into a unique expanded 
jammed state. Expansion is limited by the same inter- 
locking of polymer loops that leads to entanglements in 
polymer melts. However, the average length of straight 



segments is only N e /3 due to the random nature of the 
network. Interchain friction prevents disentanglement of 
the loops. As in dense jammed systems, the distribution 
of forces is exponential. However, the forces are tensile 
rather than compressive. We hope that contrasting this 
new jammed state with conventional ones will help un- 
ravel the microscopic origin of their common features. 
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